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Abstract. We compare the observational and theoretical spectra of the Av — 2 CO bands in a range of M dwarfs. 
We investigate the dependence of theoretical spectra on effective temperatures as well as carbon abundance. In 
general we find that the synthetic CO bands fit the observed data extremely well and are excellent diagnostics. In 
particular the synthetic spectra reasonably match observations and the best fit temperatures are similar to those 
found by empirical methods. We also examine the 12 C/ 13 C isotopic ratio. We find that fundamental 13 CO bands 
around 2.345 and 2.375 fim are good discriminators for the 12 C/ 13 C ratio in M dwarfs. The 2.375 /an is more 
useful because it doesn't suffer such serious contamination by water vapour transitions. Our current dataset does 
not quite have the wavelength coverage to perform a reliable determination of the 12 C/ 13 C ratio in M dwarfs. 
For this we recommend observing the region 2.31-2.40 /im at a resolution of better than 1000. Alternatively the 
observational problems of contamination by water vapour at 2.345 /im maybe solved by observing at resolutions 
of around 50000. We also investigated the possibility of using the Av = 1 CO bands around 4.5 fim. We find that 
the contamination due to water vapour is even more of a problem at these wavelengths. 
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1. Introduction vibrational levels of two-nuclear molecules CO is much 



More than 70% of stars in the vicinity of the Sun are M 
dwarfs. These numerous low-mass stars (0.08 M Q < M < 
0.6 M Q ), together with substellar objects - brown dwarfs 
(M < 0.08 M Q ) can contain an appreciable amount of 
the baryonic matter in the Galaxy. Research of M dwarf 
spectra are of interest for many branches of modern astro- 
physics. Verification of the theory of stellar evolution and 
structure of stars, the detection among M dwarfs of a sub- 
set of young brown dwarfs, the physical state of plasma 
of their atmospheres at low temperatures, as well as the 
chemical and physical processes of dust formation are only 
a few of them. 




The dominant opacity sources in the optical and IR 
spectra of M dwarfs is absorption by band systems of 
diatomic molecules, such as TiO and VO, as well as 
rotational- vibrational bands of H2O . M dwarf infrared 
spectra additionally contain absorption bands of CO. 
Apart from being readily identifiable in M dwarf spec- 
tra, CO bands have the advantage of being extremely 
well modelled. This is in contrast to the relatively poor 
modelling of other molecules and atomic lines in M dwarf 
atmospheres. For example, the structure of rotational- 
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12 C/ 13 C ratio in M dwarf atmospheres is especially inter- 
esting. Following the conventional theory of stellar evolu- 
tion (see Aller & McLaughlin 1965), M dwarfs save their 
initial 12 C/ 13 C from their time of formation. Since the 
galactic 12 C/ 13 C ratio is expected to change by around 
a factor of four over the lifetime of our Galaxy, the de- 
termination of the 12 C/ 13 C ratio for M dwarfs potentially 
gives a strong constraint on their age. However, deter- 
mining the 12C/ 13C ratio is only a strong constraint on 
age if that ratio is a single-valued function of time and 
Galactic location. If that underlying assumption is not 
valid, any dispersion of 12 C/ 13 C in M-dwarf atmospheres 
might give for us some clues about mixing processes inside 
our Galaxy. Until now, measurements of the 12 C/ 13 C ra- 
tio for planetary nebulae, red giants and asymptotic gi- 
ant branch stars do not fit particularly well with theory. 
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Nonetheless, the discrepancies found for such evolved stars 
may well arise from the lack of realistic stellar models. 
Atmospheric models need to include additional physical 
processes in their prescription for mixing between nucle- 
osynthetic cores and observable atmospheres (Palla et al. 
2000). M dwarfs, on the other hand, are not expected to 
modify their 12 C/ 13 C and are fully convective. Thus the 
12 C/ 13 C ratios for a diverse sample of M dwarfs is ex- 
pected to be a relatively straight-forward function of time 
and galactic location. 

Until now studies of CO bands in the K-band were 
mainly used for the quantitative analysis of the evolution 
of 12 C/ 13 C in atmospheres of red giants (cf. Hinkle 1978; 
Sneden & Pilachowski 1984; McGregor 1987; Lazaro et al. 
1991; Suntzcff & Smith 1992; Scott et al. 1994; Pilachovski 
et al. 1997). However, they are also observed in the spec- 
tra of many objects - from symbiotic stars (Schild et al. 
1992), Be stars (Kraus et al. 2000) and supernova rem- 
nants (Gerardy et al. 2000) up to planets (Forbes et al. 
1970). Baldwin et al. (1973) and many others since have 
shown the prevalence of the CO bands in the K band 
spectra of M dwarfs. Viti et al. (2002) have shown that 
the CO bands in the K band are relatively well-matched 
by synthetic spectra and thus can be used to determine 
M dwarf properties. In our study we make comparisons 
of synthetic and observed spectra across the 2.3-2.4 /j,m 
region to quantify the usefulness of the CO bands as a 
diagnostic of effective temperatures and metallicities for 
M dwarfs. 



2. Observations 

The properties of observed stars for which spectra are fit- 
ted in our paper are listed in Table 1. Observations were 
made with the Cooled Grating Spectrometer 4 on the UK 
Infrared Telescope (UKIRT) on Mauna Kea, Hawaii. The 
instrument then had a 58 x 62 InSb array which was 
moved in the focal plane in order to three times over- 
sample the spectrum. Sky subtraction was performed by 
nodding the telescope approximately 30 arcsec up and 
down the slit, ensuring that during alternate 'object' and 
'sky' observations the star remained on the detector. The 
observations presented in this paper were made during 
three nights, 1992 May 8, 10 and 12 in reasonable opti- 
cal seeing (0.75-1.5 arcsec) and of atmospheric humidity 
(10-50 per cent). The airmass difference between object 
and standard never exceeded 0.05 and so we are confident 
that the spectra have good cancellation of atmospheric 
features. 

The 150 lines mm -1 grating was used in third order 
with the 150 mm focal length camera at a central grating 
wavelength of 2.34 /im. This grating position was chosen 
to be closed to the centre of the CO bands, at a maximum 
in the grating efficiency have high atmospheric transmis- 
sion (>95 per cent) and gives a relatively high resolution, 
A/AA, where A A is the detector resolution at wavelength 
A, of 1085 (equivalent to 276 km/s). 



2.1. Standards 

Stars in the spectral type ranges B5-A5 and F6-G0 were 
used to remove the effects of atmospheric absorption. 
These standards were expected to be featureless at the 
spectral resolution used and to be well described by a 
Raylcigh-Jeans tail. However, the F6-G0 standards were 
not very useful as they show weak metal lines. Although 
the lines tended to be weaker and different from the M 
dwarf features we avoided using them where possible. In 
the reduction of spectra, we have followed the same pro- 
cedures as in Jones et al. (1996). All observations were 
wavelength calibrated using lines from observations of an 
argon lamp in the CGS4 calibration unit. This procedure 
is typically accurate to 0.1 AA. 

3. Procedure 

We used local thermal equilibrium (LTE) model atmo- 
spheres of M dwarfs with effective temperatures T c ff = 
2400-3800 K from the NextGen grid of Hauschildt et al. 
(1999). Unless otherwise mentioned all models are for log 
g = 5.0 and solar metallicity. Computations of LTE syn- 
thetic spectra were carried out by the program WITA6 
(Pavlenko 2000) assuming LTE, hydrostatic equilibrium 
for an one-dimensional model atmosphere and without 
sources and sinks of energy. The equations of ionisation- 
dissociation equilibrium were solved for media consisting 
of atoms, ions and molecules. We took into account ~ 100 
components (Pavlenko 2000). The constants for equations 
of chemical balance were taken from Tsuji (1973). We 
adopt the dissociation potential D =9.5 eV for H 2 and 
D =11.105 cV for 12 C 16 , 13 C 16 . It is worth not- 
ing that the chemical balance in M dwarf atmospheres is 
governed by the CO molecule (see section 4.1). 

Line lists of 1 H), 6 were computed using the AMES 
database (Partrige & Schwenke 1998). The partition func- 
tions of H 2 were also computed on these data (Pavlenko 
2002). We used 12 C 16 and 13 C 16 line lists of Goorvitch 
(1994). The CO partition functions were taken from 
Gurvitz et al. (1989). The atomic line list was taken from 
VALD (Kupka et al. 1999). 

The profiles of molecular and atomic lines are deter- 
mined using the Voigt function H(a,v), parameters of 
their natural broadening C'2 and van der Waals broad- 
ening C4 from databases (Kupka et al. 1999) or in their 
absence computed following Unsold (1955). Owing to the 
low temperatures in M dwarf atmospheres and conse- 
quently, electron densities, Stark broadening may be ne- 
glected. As a whole the effects of pressure broadening pre- 
vail. Computations for synthetic spectra were carried out 
with a step 0.00005 /<m for microturbulent velocity v t = 
2 km/s. 

The instrumental broadening was modelled by trian- 
gular profiles set to the resolution of the observed spec- 
tra. To find the best fits to observed spectra we fol- 
low the scheme of Jones et al. (2002). Namely, for every 
spectrum we carry out the minimisation of a 3D func- 
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Table 1. The properties of observed stars are taken from Jones et. al. (1996). Spectral types are taken from Kirkpatrick, 
Henry & McCarthy (1991) 



Object 


d (pc) 


KIN 


COL 


Sp Type 


Mk(CIT) 


Vc-Kctt 


Iccd-L' mko 


Gl 411 


2.5 


OD 


OD/H 


dM2 


6.34 


4.11 


2.24 


Gl 699 


2.39 


OD/H 


OD/H 


dM4 


8.20 


5.04 


2.59 


Gl 406 


1.83 


OD 




dM6 


9.19 


7.37 


3.70 


VB 10 


5.79 


OD 




dM8 


9.99 


8.70 


4.65 


LHS 2924 


10.72 


OD 




dM9 


10.52 


8.91 


5.19 


tion S = 


f(Xs,Xf 


%w ) — 


1/N x 


E(l - Fobs 


/ F syn t) 2 , 


20 r 









where F b s ,F syn t are observed and computed fluxes, N is 
the number of points in observed spectrum to be fitted, 
relative shift in wavelength scale, a normal- 
isation factor which was used to coincide observed and 
computed spectra and parameter of instrumental broad- 
ening, respectively. 

4. Results 

4.1. CO and H2O molecules in M dwarf atmospheres 

The overall chemical balance in the atmospheres of cool 
stars depends significantly on the C/O = log N(C) - log 
N(0)(see Tsuji 1973). In the atmospheres of cool (T ef f < 
3600 K) oxygen rich stars C/O < 1, almost all atoms of 
carbon are bound in CO molecules. CO and H2O are the 
most abundant molecules containing oxygen atoms. 

In the atmospheres of M dwarfs with T ff > 2600 K, 
molecules of CO dominate by number in comparison 
with other oxygen containing molecules (Tsuji 1973). 
Throughout most of the atmosphere changes of the ra- 
tio C/O by ± 0.2 dex have little effect on the molecular 
densities of CO (Fig. |l|). In contrast the dependence of 
molecular densities on H2O on log N (C) is more com- 
plex. At lower T c ff (e.g. Fig. lb) the dependence of nn 2 
on log N(C) in M dwarfs atmospheres weakens apprecia- 
bly though the chemical balance remains governed by CO. 
In general, H2O and CO are competitors for bound-free 
oxygen atoms. Therefore their molecular densities respond 
in opposite ways to changes of carbon abundance. For so- 
lar abundances of oxygen and carbon molecular densities 
of H2O and CO are comparable. CO molecular densities 
increase only a little when log N(C) increases from -3.48 
up to -3.12^|. At the same time CO reduces molecular den- 
sities of H2O . However, relative changes of H2O densities 
are comparatively weak for variations of log N(C) in the 
range from -3.48 to -3.18. Only when the carbon abun- 
dance rises to log N(C) ~ -3.18, molecular densities of 
H 2 decrease, as an appreciable fraction of oxygen atoms 
become bound in CO. 

The relative strengths of H2O and CO bands at 2600 
K and 3600 K are shown in Fig.||. These spectra were 
computed separately by switching on/off the absorption 
of H 2 or CO bands. 
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Fig. 1. Molecular densities of H2O and CO in atmo- 
spheres of M dwarfs 3600/5.0 and 2600/5.0 as a function 
of gas pressure P g . The case of log N(C) = -3.48 and log 
N(O) = -3.12 represent the solar abundances case. 



The relative intensity of the CO bands decrease from 
Aw = 1 to Av = 3. Although the second overtone bands 
(Aw = 2) of CO have intermediate intensities they are 
located in a relatively accessible region for observation. 

4.1.1. Bands of the second overtone of CO from 
2.28-2.4 //m 



In our paper we use an abundance scale EiVi = 1. 



The relative intensities of 12 C 16 , 13 C 16 and H 2 bands 
from 2.28-2.38 /im are shown in Fig. ||. These computa- 
tions are carried out for a model atmosphere with 2600 K 
and 12 C/ 13 C = 0.5. The large ratio of 12 C/ 13 C was cho- 
sen to show the positions of the main details which might 
be observed in M dwarf spectra. Overall, H 2 absorption 
dominates the infrared spectral region. However, strong 
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Fig. 2. The relative strength of 12 C 16 and H 2 bands 
in a 2600 K synthetic spectrum. The positions and relative 
strengths of the CO bands can be seen. 




2.32 2.34 2.36 

Wavelength in microns 



Fig. 3. The relative strength of 12 C 16 , H 2 and 
13 C 16 bands in the spectrum of VB10 (top), and a com- 
bined spectra computed for the case of 12 C 16 + H2O and 
12 C 16 Q + i3 C i6 + jj 2 q absorptions. The model atmo- 
spheres shown are 2600 K, 12 C/ 13 C = 0.5 (bottom). The 
positions of 13 C 16 bands heads are shown by arrows. 



features are also formed there by first overtone bands 
of 12 C 16 . The positions of the 13 C 16 bandheads are 
marked by vertical arrows in Fig. 3. The results of our 
13 C 16 modelling of spectra of our dwarfs are discussed 
in sections 4.1.3 and 4.1.4. 
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log N(C) = -3.18 
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LHS 2924 
2600/5.0, log N(C)=-3.48 
2800/5.0, log N(C)=-3.28 




2.31 2.32 2.33 2.34 2.35 2.36 2.37 
Wavelength in microns 



Fig. 4. Top: The dependence of S on T cff and log N(C). 
Bottom: 2600 and 2800 K synthetic spectra compared to 
the observed LHS2924 spectrum. 



4.1.2. Fits to the observed spectra 

Our numerical procedure for the determination of best fits 
allows us to quantify the possible solutions. For each ob- 
served spectra we compute synthetic spectra for different 
T e ff . We also allow for small variations of carbon abun- 
dances. Based on Fig. 1, we suggest that small variations 
of log N(C) from -3.48 to -3.28 cannot substantially af- 
fect the temperature structure of the model atmospheres. 
From the family of curves S = f(T e s > logN (C)) we choose 
"the best solution" for each observed spectrum. 

LHS2924. The observed spectrum of this M9V ob- 
ject is of relatively low signal-to-noise. For example, our 
best fit minimisation value S m i n for LHS2924 is a factor 
10 larger than VB10 (fig. ||). Our best fit provides T cff = 
2800 K for log N(C) = -3.28 and T cff = 2600 K for a solar 
abundance of carbon (log N(C)=-3.48). 

VB10. VB10 is an archetypal late-type M dwarf with 
an effective temperature T c g ~ 2600 K (Jones et al. 
1996). Lithium lines are absent from VBlO's spectrum 
(Schweitzer et al. 1996; Martin 1999) and it has typical 
old disk properties (Leggett 1992) and hence VB10 is a 
star right at the end of the main sequence. For the solar 
abundance of carbon and other elements S m i n occurs for 
2600 K. However, if we increase log N(C), S decreases, 
and we find a best fit value of 2800 K for log N(C) = 
-3.28±0.1. In general, for higher log N(C), CO bands in- 
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2800/5.0, log N(C)— 3.48 
3000/5.0, log N(C)=-3.28 




2.32 2.34 
Wavelength in microns 



Fig. 5. Top: The dependence of S on T off and log N(C). 
Bottom: The best fit of synthetic spectra to the observed 
VB10 spectrum. 



tensities increase and best fits move to higher T e g . The 
same effect is found for all dwarfs of our sample. 

GJ406. Best fits are found for 2800 K for the solar 
abundance case, and 3000 K for log N(C) = -3.28 (Fig. g). 
Our new estimation of effective temperature corresponds 
better with empirical values for the effective temperature 
of the spectral class M6V than the Jones et al. (2002) 
analysis found using a similar technique but using a re- 
gion dominated by water vapour. Our best fits for GJ406 
are found for solar metallicity rather than the metal poor 
result found by Jones et al. (2002). 

GJ699. The observed spectrum of GJ699 (M4V) is of 
lower resolution and taken from the setup used in Jones et 
al. (1994). The details in observed spectrum are best fit by 
metal-poor (log N(C) = -3.68; i.e., [C] = -0.2) synthetic 
spectra of T ff = 3200 K (see fig. 0). The same best fit 
temperature is also found for the solar abundances log 
N(C)= -3.48, ([C] =0) case. However, a carbon-deficient 
case log N(C) = -3.68 QC]=-0.2) provides lower S. We note 
that Jones et al. (2002) found an effective temperature 
T off = 3300 K and metallicity [M/H] = - 0.5 ([O] = - 
0.5). On the other hand, our fits with log N(C) = -3.88 
([C] = -0.4) provide less satisfactory fits, i.e fits of higher 
S. We have probably reached the limit at which a single 
metallicity parameter is appropriate for fits to the spectra 
of M- dwarfs. 

GJ411 is the hottest dwarf in our sample. The fit 
to the observed spectrum is shown in Fig. 0. Within our 



Fig. 6. The same as in Fig. | but for GJ406. Synthetic 
spectra are shown for T e a = 2800 and 3000 K. 



log N(C)=-3.88 
log N(C) = -3.68 
log N(C) = -3.48 
log N(C) = -3.28 
log N(C) = -3.18 



3200 



GJ 699 
3200/5.0, log N(C)=-3.48 
3200/5.0, log N(C)=-3.68 
3300/5.0/-0.5 




2.31 2.32 2.33 2.34 2.35 
Wavelength in microns 



2.36 2.37 



Fig. 7. The same as in Fig. | but for GJ699. Synthetic 
spectra computed for 3400 and 3600 K synthetic spectra 
are shown. 
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log N(C) = -3.68 
log N(C) = -3.48 
log N(C) = -3.28 
log N(C) = -3.18 




GJ 41 1 
3400/5.0, log N(C)=- 3.48 
3600/5.0, log N(C)=-3.68 



3700 3800 





2.342 2.3425 2.343 2.3435 2.344 2.3445 2.345 2.3455 2.346 2.3465 2.347 
Wavelength in microns 




2.29 2.3 2.31 2.32 2.33 2.34 2.35 2.36 2.37 
Wavelength in microns 



2.34 2.345 2.35 2.355 2.36 2.365 2.37 2.375 2.38 2.385 
Wavelength in microns 



Fig. 8. The same as in Fig. 0, but for GJ411 including 
synthetic spectra for 3400 and 3600 K. 

observed spectral region the intensities of observed and 
computed band heads of are fitted for T ff = 3400 K 
and log N(C) = -3.48 (Fig. |). Although the relatively 
poorer description of the M dwarf spectrum slope testi- 
fies to the presence of certain problems in the definition 
of its physical characteristics: abundances, homogeneity of 
atmospheres, chromospheres or other reasons (see Section 
5). 

4.1.3. Isotopic ratio 12 C/ 13 C in M dwarf atmospheres 

The temperatures in M dwarf interiors are relatively low 
(T < 10 6 K), so that hydrogen burning is only possible 
in pp chain reactions and the hydrogen burning lifetimes 
are rather long (e.g. Laughlin et al. 1997). Thus despite 
the large convective envelopes extending down to the hy- 
drogen burning region in a short timescales, carbon iso- 
tope abundances cannot be changed even in Hubble-scale 
times. M dwarfs are expected to preserve their initial car- 
bon isotopes abundance from their time of formation. The 
situation is very different for massive (M > 0.8 M Q ) stars, 
which burn hydrogen in CNO reactions and are observed 
to change their isotopic abundances of C,N,0. The deter- 
mination of the 12 C/ 13 C ratio in the atmospheres of M 
dwarfs might provide important constraints on the theory 
of the evolution of the Galaxy. 

We investigated the determination of the carbon iso- 
topic ratio in the atmospheres of our observed M dwarfs. 
The analysis of the weak 13 C 16 bands in our compar- 



Fig. 9. Top: Spectra computed for model atmosphere 
2800/5.0, log N(C)=-3.28, for 12 C/ 13 C = 3 to 90 at a reso- 
lution of around 46000. Bottom: Fits of computed spectra 
with different 12 C/ 13 C to VB10 observed spectrum. 

atively low-resolution spectra of dwarfs is rather difficult 
(e.g., Fig. § for VB10). The CO band head at 2.345 is 
mimiced by H2O absorption. Moreover, the overall shape 
of the energy distribution at these wavelengths is governed 
by H2O absorption, Fig. |^. From the fit of our spectra for 
VB10 we can estimate the lower limit 12 C/ 13 C > 15. For 
more accurate determination using this feature it will be 
necessary to use spectra of better resolution (> 50000 for 
12 C/ 13 C in range 40 — 90, e.g., upper panel of Fig. 
The other band of 13 C 16 at 2.375 /zm provides a more 
sensitive tool for 12 C/ 13 C determination in M-dwarf spec- 
tra. From the fit of the observed band to our models we 
find 12 C/ 13 C ~ 10. Unfortunately our confidence in this fit 
is weakened because our data do not extend beyond 2.38 
/im. To be sure we are measuring this feature we need to 
have data beyond 2.38 /im. For the other M dwarfs the 
fits are not so good so it is necessary to confirm this re- 
sult with spectra of higher resolution and signal-to-noise. 
Nevertheless, from the fit of our spectra we can estimate 
the lower limit 12 C/ 13 C > 10 or all stars. 

4.1.4. Bands of the first overtone of CO from 4.3 to 
4.6 microns 

As with the Au — 2 CO bands at 2.3/im it might be 
feasible to pursue the 12 C/ 13 C ratio using the Az/=1 CO 
bands at 4.4/im. A theoretical picture of formation of the 
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Table 2. T c ff of our sample for solar metallicity compared with other methods: water vapour (Jones et al. 2002), 
infrared flux method (Tsuji et al. 1997), constant opacity (Jones et al. 1994), model fit based on temperatures derived 
by Leggett et al. (2000, 2001). The temperatures in italics are those derived for stars of similar spectral type. 



Object 


Spectral 
Type 


T c ff 
This work 


T c ff 
Water vapour 


T c ff 
IRFM 


T c ff 

Constant opacity 


T ff 
Model fit 


GJ 411 


M2V 


3400 




3510 


3471 


3450 


GJ 699 


M4V 


3200 


3300 


3210 


3095 


3100 


GJ 406 


M6V 


2800 


3000 


2800 


2670 


2600 


VB10 


M8 


2600 




2250 


2506 


2250 


LHS2924 


M9 


2600 




2120 


2219 


2125 




4.25 4.3 4.35 4.4 4.45 4.5 4.55 4.6 

Wavelength in microns 




4.36 4.38 4.4 4.42 4.44 

Wavelength in microns 

Fig. 10. Top: relative strength of 12 C 16 , H 2 and 
13 C 16 bands in theoretical spectrum of M-dwarf 
2800/5.0/0. and a combined spectra computed for the case 
of 12 C 16 + H 2 and 12 C 16 + 13 C 16 + H 2 absorp- 
tion. Bottom: synthetical spectra computed for 2800/5.0 
model atmospheres and different 12 C/ 13 C . The resolution 
R = 4600 was adopted. 



absorption features in the region of the bands of the first 
overtone of CO is shown in Fig. [l(]. Again, as for the 13 CO 
feature at 2.345 /im, lines of H2O dominate in region (see 
also Fig. |). Bands of 12 CO (and 13 CO) are actually of 
lower intensities relative to water. The band heads of 13 C 
fall such that their heads lie in the tails of 12 C bands, 
which are stronger. Therefore the determination of the 
12 C/ 13 C ratio for At/ = 1 CO bands is also problematic. 



5. Discussion 

Our effective temperature determinations are based pri- 
marily on 12 C 16 lines but also affected by H 2 lines. In 
Table 2 it can be seen that we find similar temperatures 
to those found for M dwarfs for H2O lines alone. These 
results support the expected quality of the 12 C 16 and 
H 2 line lists. 

In general, we find T e fj of our dwarfs with internal 
errors of around ± 150 K. Table 2 shows T e ff derived 
for M dwarfs by different methods. It illustrates indicates 
that whilst there is some agreement for the earlier type M 
dwarfs, the temperatures we derive for late spectral types 
are much too hot. This discrepancy arises because in the 
atmospheres of red dwarfs of spectral class M6 and later 
(Tsuji, Ohnaka & Aoki 1996; Jones & Tsuji 1997) appre- 
ciable amounts of dust forms (e.g., Tsuji 2002) that can 
absorb and/or scatter radiation. This dust changes the 
structure of model atmospheres of M dwarfs, and, hence, 
their spectra. The outermost layers of late M dwarf atmo- 
spheres increase due to appearance there of the "dusty" 
opacity. In general, we get the much more sophisticated 
picture of formation of late dwarf spectra in the "dusty" 
atmospheres in comparison with "non-dusty" case. For the 
case of our fits the dust in cooler stars causes the CO 
bands to be relatively weaker and thus best-fit by a high 
temperature, non-dusty model. Analysis of our data to 
determine reliable temperatures for VB 10 and LHS 2924 
thus requires the use of dusty models. Such models are 
now in preparation (Pavlenko et al., in preparation). 

NLTE (Non Local Thermal Equilibrium) effects pro- 
vide another possible source of error. Carbon et al. (1976) 
investigated NLTE effects in CO lines of the fundamental 
systems formed in the atmospheres of stars of giant late- 
type spectral classes. However, they used a rather simple 
model of rotational- vibrational levels of CO. Furthermore 
they considered the model atmospheres of red giants, 
where the densities are much lower, than in M dwarf at- 
mospheres. Nevertheless, they showed the possibility of 
NLTE effects in CO lines, formed on depths ri Mm < 10~ 4 . 
CO lines in our cases are strong, they are formed deep 
in the atmosphere (due to the low opacities in contin- 
uum). The formation of absorption lines of CO occurs in 
M dwarfs atmospheres on a background of strong H 2 O ab- 



8 



Yakiv V. Pavlenko and Hugh R.A. Jones: CO bands in M dwarf spectra 



sorption. The additional absorption should strengthen the 
thermalisation of the radiation field on CO lines. The tem- 
perature structure of M dwarf atmospheres is smoother, 
in comparison to the giants. Indeed, their structure is 
governed by convection, which "smoothes" the temper- 
ature contrast between photosphere and upper layers, 
which are more pronounced in atmospheres of red gi- 
ants. Atmospheres of cool dwarfs are denser, than giants 
which strengthens the thermalisation of the radiation field. 
All these factors should reduce NLTE effects in CO lines 
formed in M dwarf atmospheres. More precise estimations 
are possible and would be desirable when making compar- 
isons with higher resolution data though we consider to be 
beyond the scope of the current work. 

Our computations were carried out for model atmo- 
spheres of solar abundance. Such abundances should be 
prevalent for the majority M dwarfs in the vicinity of solar 
system. Although, our sample of M dwarfs is rather small, 
we have found some deviation in carbon abundances from 
the solar case. A more refined solution should be made 
using self-consistent abundances of oxygen and nitrogen. 

The model atmospheres used in this work were com- 
puted in the framework of the classical assumptions by 
Hauschildt ct al. (1999). Energy transfer in atmospheres 
is provided by radiation field and convection. The relia- 
bility of the computed model atmosphere is determined 
by both the completeness of the opacity sources as well 
as the model of convective processes(Ludwig et al. 2002). 
Absorption by H 2 bands determines the opacity in the 
infrared spectra of red stars. However, the existing line 
lists are not sufficiently complete (e.g., Jones et al. 2002). 
The CO line list used in this work was computed for so- 
lar work and so is believed to be more than adequate for 
the lower energy states accessed in M dwarf atmospheres 
(Goorvitch 1994). Overall we do not consider that incom- 
pleteness in the H 2 line lists will cause temperatures 
derived from CO bands to change by more than 150 K. 

The spectra presented here can not be used to de- 
termine reliable 12 C/ 13 C ratios. Nonetheless we propose 
that with a slightly different wavelength coverage than ob- 
served for this study, the 13 CO feature at 2.375 /xm should 
be a strong constraint on the 12 C/ 13 C ratio. For the other 
13 CO feature we investigated the situation is complicated 
by 13 C 16 bands forming amongst strong H 2 absorption. 
One way of dealing with such 'contamination' would be to 
obtain higher resolution data. In the upper panel of Fig. 9 
we show synthetic spectra at a resolution of around 46000, 
here even at low 12 C/ 13 C ratios the 12 CO and 13 CO bands 
can be distinguished. We have also investigated the possi- 
ble use of the 4.4-4.6 micron region. However, this region 
does not appear to improve the possibility of determining 
12 C/ 13 C . 
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